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ABSTRACT 

The interpretation of supernova (SN) spectra is essential for deriving SN ejecta proper- 
ties such as density and composition, which in turn can teh us about their progenitors 
and the explosion mechanism. A very large number of atomic processes are important 
for spectrum formation. Several tools for calculating SN spectra exist, but they mainly 
focus on the very early or late epochs. The intermediate phase, which requires a NLTE 
treatment of radiation transport has rarely been studied. 

In this paper we present a new SN radiation transport code, nero, which can look 
at those epochs. All the atomic processes are treated in full NLTE, under a steady- 
state assumption. This is a valid approach between roughly 50 and 500 days after the 
explosion depending on SN type. This covers the post-maximum photospheric and the 
early and the intermediate nebular phase. 

As a test, we compare NERO to the radiation transp ort code of I Jerkstrand et all 

(|2011h and to the nebular code of iMazzali et al.l ()2001l) . All three codes have been 
developed independently and a comparison provides a valuable opportunity to inves- 
tigate their reliability. Currently, nero is one-dimensional and can be used for pre- 
dicting spectra of synthetic explosion models or for deriving SN properties by spectral 
modelling. To demonstrate this, we study the spectra of the 'normal' SN la 2005cf 
between 50 and 350 days after the explosion and identify most of the common SN la 
line features at post maximum epochs. 
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1 INTRODUCTION 

Supernovae are classified by their spectra. T he most com- 
mon types are SNe la, lb, Ic, lib and II (e.g. iTuratto et al.l 
I2OO7I ). This classification links directly to the explosion 
mechanism, which can be either a thermonuclear explosion 
or the collapse of a massive star. The connection can be 
established since the progenitor and the explosion mecha- 
nism leave a unique imprint on the resulting SNe, which is 
refiected by their spectra. 

In turn, one can learn about the properties of observed 
SNe by spectral analysis. The earliest attempts to interpret 
the spectra of SNe reach back to the very beginning of SN 
astronomy, when SNe I were classified as hydrogen poor and 



SNe II as hydrogen rich (e.g. lMinkowskilll94lh . Since then, 
the interpretation of SN spectra has been steadily refined 
and several advanced tools to study SN spectra formation 
exist today (see below). 

SN spectra are usually grouped into photospheric and 
nebular. It is not trivial to distinguish between the two 
phases on physical grounds since most radiation process 
important during the earliest phases play also a role at 
late times. However, nebular phase spectra are dominated 
by clear forbidden-line emission features, which become 
dominant at about 200 days after the explosion. Although 
photo- excitation a nd -ionisation proce s ses are important at 
later epochs (e.g. Ibi fc McCravl Il996l : Ide Kool et all Il998l : 
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Ijerkstrand et al1l201ll . also see this paper), they are abso- 
lutely dominant at earlier epochs. 

Directly after the explosion the gas is in local ther- 
modynamic equilibrium (LTE), which simplifies radiation 
transport. However, shortly after maximum light at latest 
non-local thermodynamic equilibrium (NLTE) effects be- 
come important. Very little attention has been paid to these 
intermediate epochs so far, which require a full NLTE treat- 
ment of radiation transport. 

Codes treating the early pha se hav e for exampl e 



been developed b y [Mazzah fc Lucy! (Il993l ): ILucvI (| 19991 ): 
iKasen et af] (120061 ) : lKromer fc SimI (|2009l ') . These early time 
codes rely at least partially on the LTE assumption and do 
not treat forbidden-line emission. For later epochs numeri 
cal tre a tments have for example been d eveloped by Axelrod 



198C) 



1992 



Ruiz-Lapuente fc Lucvl (1 19921): iKozrna fc Fransson 



1998h: lEastman fc Pintol lll993h: 



1998h : lMazzah et all (|200ll ): |jerkstrand et al 



de Kool et al 
(|201ll ).Inad 



dition to thos e spectral codes specified for SNe, there ex- 
ist others (e.g.lPauldrach et al.lll99^ : iHauschildt et al.ll 19971 : 
iHiUier fc Milleij|l998l ) which can be used for calculating SN 
spectra. The list presented here is a selection and is far from 
complete. 

Here we present a new non-thermal equilibrium radia- 
tion transport (nero) code which addresses especially the 
intermediate epochs at about 50 — 200 days after the explo- 
sion. We can also treat later epochs, as long as the gas shows 
no strong deviations from steady-state (see Section [SJ. 

In Section [2] we describe the new code. In Section [3] 
we co mpare NERO to th e stead y-state radiation transport 
code of Jerkstran d et al] |20l3) and to the nebular code of 
iMazzali e t al. (200"^ We also compare synth etic spectra ob- 
taine d by using NERO on a SN la W7 model l|Nomoto et al.l 
1 1984 ) with observations of the proto-typical 'normal' SN la 
2005cf at epochs between 50 and 350 days after the explo- 
sion. In Section |4] our results are discussed. 



2 CODE DESCRIPTION 

NERO is built on the assumption that the SN ejecta are in 
a steady-state, i.e. heating and ionisation are balanced by 
cooling and recombination, which is valid between roughly 
50 — 500 days after the explosion (e.g. lAxelrodlll980l ). These 
limits however depend strongly on the ejecta density and 
composition. Under this assumption the spectral emission 
can be calculated at any epoch without considering the radi- 
ation history of the gas. The phase we can treat begins when 
the SN light-curve starts to follow the decay of Ni and Co, 
which means that the gas is in thermal equilibrium and ends 
when the gas falls out of ionisation equil ibrium or when adi- 
abatic cooling becomes important (e.g. iKozma fc FranssonI 
Il992l ). 

In a first step, the energy deposition from the decay of 
radioactives is calculated, from which we derive the non- 
thermal electron ionisation and excitation rates. Using an 
initial guess for the ionisation and excitation state and the 
temperature of the gas we derive its new ionisation and ex- 
citation state by solving the statistical equilibrium equa- 
tions. From the ion fractions and occupation numbers one 
can calculate the radiation field. This radiation field is trans- 
ported through the SN ejecta using a ray-tracing transport 
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Figure 1. Code scheme of NERO. 



scheme, from which we derive photo-ionisation and excita- 
tion rates. The ionisation and excitation states are iterated 
until the electron density and temperature and the radiation 
field have converged. A scheme of NERO is shown in Figure 

m 

Currently, NERO is one-dimensional, but a three- 
dimensional version may be available in the future. For a 
one-dimensional model consisting of about 20 radial shells, 
one calculation takes a few minutes on a standard desktop 
computer, depending on the composition and density of the 
model and the desired resolution. All synthetic spectra ob- 
tained for this paper using NERO required a calculation time 
of the order of 5 minutes, which is excellent for a fully-fledged 
NLTE radiation transport code. Still, we plan to parallelise 
and further optimise NERO in the future, to make modelling 
of observed SN spectra as efhcient as possible. 



2.1 Radioactive deposition 

NERO can treat radioactive deposition by ^^'^'^Co and by 
The half-lives of these is otopes can for example be 
found in ISeitenzahl et al.l l|2009l ). Including other radioac- 
tives would be simple. The energy deposition process is 
treated with a Monte Carlo approach, where Compton scat- 
tering, pair-creation and photo-ionis ation are taken i nto ac - 
count, as for example described in ISim fc Mazzalil (|2008l ). 
We calculate the spectrum of the emerging 7-radiation. 
However, there are no observations available to compare 
with. The accuracy of our deposition routine has been veri- 
fied by comparing the energy depositio n rates to calcu l ations 
perf ormed with the spectra l codes of iMazzali et al.l (1200 ll ) 
and Ijerkstrand et all l|201lf ) . From the deposited luminos- 
ity we calculat e non-thermal ionisation and e xcitation rates 
as described in Im aurer fc Mazz ali (20l3) and lMaurer et al] 
(^oioi). To evaluate the non-thermal ionisation rates we es- 
timate the fraction of the deposited luminosity causing ion- 
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isation using the Bethe approximation (e.g. lAxelrod|[l980l ). 
The non-thermal excitation rates are calculated using the 
optic al approximation (e.g. lAxelrod ll98C ; iRozsnvai et al.l 



ligsd ). It was found bv lMaurer et al.l (|2010af ) that this treat- 
ment of non-thermal ionisation and excitation is accurate to 
at least 20% for hydrogen and helium. We also compared 
our non-thermal electron ionisation and excitation rates for 
various ions to rates calculated with the radiation trans- 
port code of Ijerkstrand et al.l ( 20 111'), which make s use o f 
the Spencer- Fano approach of Kozma fc Franss"onl (|l992l ). 
We find overall good agreement, with deviations usually of 
the order of 10% and ~ 50% in the worst case. Especially for 
the iron-group elements the non-thermal rates are strongly 
affected by inaccurate or poorly known atomic data and a 
high degree of uncertainty has to be accepted. This may be 
the largest source of uncertainty in all spectral calculations 
of SNe la after maximum light. 

2.2 Ionisation &: Recombination 

The ionisation state of the gas is calculated under the as- 
sumption of statistical equilibrium, balancing non-thermal 
electron ionisation, photo-ionisation and radiative and di- 
electroni c recombinatio n with the charge-exchange reactions 
listed in ISwart^ ijlQQA . The non-thermal electron rates are 
obtained as described in Section 12.11 The photo-ionisation 
rates are obtained from ray-tracing radiation transport (de- 
scribed below) using the ground state photo-ionisation cross- 
sections from TIPTOP Bas«El and a simple approximation 
for the lowest 40 excited states. Above the correspond- 
ing ionisation thresholds, the excited-state photo-ionisation 
cross-sections are assumed to be constant fractions of the 
ground-state cross-sections. These fractions decrease with 
increasing main quantum number of the excited states. All 
rates are corrected for stimulated recombination. We do not 
treat photo-ionisation from states higher than 40 since they 
seem to have no noticeable influence on the spectra at the 
epochs of interest (50 — 500 days). We use t he radiative and 
di-elec tronic total recombination rates of iMazzotta et al.l 
lll998h. Ground state recombi nation rates are taken from 
lAldrovandi fc Pequigno"^ l| 19731) or are set to 10% of the total 
radiative recombination rate, if not available. All the elec- 
trons recombining into excited states are distributed equally 
to the lowest 40 excited states from where they undergo the 
complete NLTE process. We plan to improve our excited 
state ionisation and recombination data base in the future. 



2.3 Excitation 



We use the line data collection of iKurucz fc Beiil l|l995h . 
which roughly contains 25.000 atomic levels and 500.000 
lines. In principle, we can treat all elements from H to Ni 
and all ions from ionisation state I — III (IV is taken into ac- 
count for the ionisation equilibrium but does not contribute 
to the radiation field). However, the atomic data are poor 
for many ions. The atomic e xcitation state s are calculated 
by solving a rate matrix (e.g. lAxelrodlll980 ). including non- 
thermal electron excitation rates (see Section I2.1|l , photo- 
ionisation and excitation rates obtained from ray-tracing 



radiation transport (see below) using the Sobolev approxi- 
mation, spontaneous radiative de-excitation, recombination 
into excited states, th ermal electron (de-) excitation, co ntin- 
uum destruction (e.g. IChugai|[l987l : iLi fc McCra^ll995l ) and 
two-pho ton emission (TPE) o f H l(2s^S) a nd He l(2s^'^S) 
[e.g. see iKaplan et all l|l972D : iDrake et al.1 l|l969l ) for TPE 
rates]. 

With the electron density and temperature obtained 
from the iteration process, we calculate coUisional (de- 
) excitation rates. CoUisional data are taken from TIP- 



(e.g. iBerrington et al.l Il982l: Haves & Nussbaumer 
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iBerrington' '1988": 'Manas et al. 


1988*; 'Scholz et al 




199C 


Callawav 1994: Mclondez et al.l 


2007: Bautista et al. 


|2009|) 



Since our coUisional database is far from complete we plan 
to regularly add and update coUisional atomic data. If 
not available, the coll i sion strengths are approximated (e.g. 
Ivan Regemorte j 1 19621 : 1 AxehodI I198G| ) . However, for several 
hundreds of lines there is coUisional data from the litera- 
ture. A serious problem at intermediate epochs, especially 
for treating SNe la, is the absence of reliable coUisional data 
for Co. 



2.4 Radiation Transport 

From the ionisation and excitation states obtained in the 
previous steps we calculate a radiation field, which is repre- 
sented by a certain amount of photon packets (typically, a 
few 100.000 per shell in total). These are sent through the 
SN envelope in random directions. On their way out they 
propagate on straight lines and encounter bound-bound and 
bound-free absorption and electron-scattering. The proba- 
bilities for line scattering are calculated in the Sobolev ap- 
proximation. While transported, the photon packets lose 
parts of their energy according to the respective optical 
depths (or change their direction after electron scattering), 
and can be absorbed to 100% if the optical depth is much 
larger than one. The absorbed photon packets are re-emitted 
in random directions in the next iteration step after tak- 
ing part in the NLTE excitation matrix calculation (Section 
I2.3|l . Therefore, all the absorbed energy undergoes the full 
NLTE process, including fluorescence, up- and down-ward 
electron collisions, photo-ionisation, recombination, contin- 
uum destruction and two-photon emission. 



3 CODE RESULTS 

3.1 Comparison to other codes 

In this secti on we compare NERO to the radiation trans- 
Jerkstrand et al.l ll201ll ) [RTJ] and to the neb- 
Mazzah et al.l (120011 ) [RTM]. RTJ is a steady- 



port code of 
ular code of 

state radiation transport code. Although developed inde- 
pendently, it is built on very similar physical assumptions 
as NERO. While NERO treats the NLTE process completely, 
RTJ does not allow up- ward excitation of photo-excited lev- 
els. However, at least at the late epochs, which have been 
chosen for the code comparison, this seems to have no ob- 
servable influence on the synthetic spectra (see below) . RTJ 
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Table 1. SN 19941. (A) RTM one-zone model (B) NERO one-zone 
model (C) NERO small-scale separation model 
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0.2 
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0.03 
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0.07 


0.2 


0.05 


0.1 


0.01 


0.002 


0.05 



has be en applied to study the very late phase of SN 1987A 
so far l|Kiaer et al.ll20ld : iJerkstrand et al.ll201lD . 

R TM is a nebular code based on th e ideas of lAxelrodl 
l|l98(tl and lRuiz-Lapuente fc Lucvl (|l992h . therefore neglect- 
ing all radiation trans port effects. It has widely been used 
in th e literature (e.g. ISilverman et al.l |2009| : iMazzali et al.l 
to study nebular spectra of all types of SNe. 

While NERO and RTJ are currently available in one- 
dimensional v ersions only, there are three-dimensional ver- 
sions of RTM l|Maeda et al.ll2002l : [Maurer et al.ll2010bl ). 

A code comparison is interesting, since all three codes 
have been developed independently and use different numer- 
ical methods, physical assumptions and partially different 
atomic data. 

For the c omparison of RTJ an d NERO, we chose the 
13C model of IWooslev et all (|l994l '). At 200 and 400 days 
after the explosion excellent agreement (see Figures [5] & [31) 
is observed. The 13C model is found to show very strong 
Ca II emission in both calculations, which can however be 
explained to be a mixing effect. In the 13C model all the Ca 
is mixed with the other elements microscopically. Since Ca 
II has a low excitation potential and large collision strengths 
it radiates strongly, if it is mixed into large amounts of hy- 
drogen, helium or oxygen. 

Another comparison of NE RO and RTJ wa s performed 
for a Type la W7 model (Nomoto et al.lll984l 'l at 94 and 
338 days after the explosion. The agreement at both epochs 
is good (see Figures |4] & [5| . However, at 94 days after the 
explosion there is a strong deviation around 4600 A and 
5900 A. These features are caused by Fe ill and Co ill, re- 
spectively. The fraction of these ions is similar in both cal- 
culations, which means that the ion abundance cannot be 
a main reason for this differences. Since both features are 
dominated by coUisional excitation from ground levels and 
since the electron temperature and the density are similar in 
both calculations, it is likely that the differing sets of atomic 
data used in NERO and RTJ are responsible for most of the 
deviation. Also, there is is an important fraction (~ 15% at 
5000 km/s, increasing with velocity) of Fe iv and Co iv in 
the NERO calculation, which is neglected by RTJ at the mo- 
ment. The recombination of these ions (and photo-ionisation 
of Fe III and Co ill) influences the radiation field and the cas- 
cading of UV radiation can hardly be followed in detail. It 
is interesting to note that the Fe ill feature produced with 
RTJ is more consistent with observed SNe la spectra, while 
this is true for the Co ill feature produced with NERO (see 
below) . 

For the comp arison of RTM and NERO we chose a 
Type la W7 model ( Nomoto et"aLll 19841 ') and the 'standard' 
Type Ic SN 19941 (|sauer et alTHooel ). The comparison of 



spectra obtained with RTM and NERO for the W7 model 



jNomoto et al.|[l984h shows reasonable agreement (see Fig- 
ure [6l). The core of the W7 model, which is observed in 
the nebular phase, consists of almost pure Fe (from ^^Ni 
decay). In NERO Fe I is photo-ionised almost completely. 
The ratio of Fe 11 to Fe ill is dominated by non-thermal 
electron ionisation at late epochs, which is also treated by 
RTM. In RTM the fraction of Fe I is always set to zero. 
Therefore, the ionisation state (and hence the electron den- 
sity) of a pure Fe core obtained by NERO and RTM is sim- 
ilar for this model. The electron temperature of the Fe 
plasma obtained by NERO is about 10% lower than that ob- 
tained with RTM, because of the presence of excitation pro- 
cesses (photo-excitation, excitation by excited-state recom- 
bination, non-thermal electron excitation, more coUisional 
transitions) which are neglected in RTM calculations. 

Since RTM has exclusively been used to derive SN core 
ejecta properties (and not for predicting spectra from ex- 
plosion models) , we also compare modelling results of RTM 
and NERO. For this compariso n, we chose SN Ic 1 9941, which 
has been studied in detail bv ISauer et al.1 l|2006l ). They de- 
rived a total core mass (v < 5500 km/s) of 0.43 M© and a 
^^Ni mass of 0.07 M©, which was also found to be consis- 
tent with the observed hght curve of SN 19941. We fit the 
observed spectr um with one - zone models (see Figure [T}, as 
it was done by ISauer et all (120061 ) [see Table [H model 'A' 
(RTM) & 'B' (nero)]. 

In general, the masses estimated for Na, Mg, Si, S 
and Ca deviate in RTM and NERO calculations. This is 
expected, since these elements are strongly influenced by 
photo-ionisation, which is not treated in RTM. Therefore, 
one derives more mass for the elements which are esti- 
mated from the neutral component (e.g. Na, Mg) and less 
for those which are estimated from ionised states (e.g. Ca) 
using NERO. 

Another important difference is found for the estimate 
of the total and the ®^Ni mass. These two quantities can be 
considered as the main properties of any SN and should be 
in unison with its hght curve. While the RTM calculat ion is 
consistent with the hght curve (see ISauer et al.ll2006l) . the 
NERO calculation seems to show too much total mass and 
too little '^'^Ni (see Table[T] model ' B') to be con sistent with 
the light curve modelling results of ISauer et al.l (2006). 

This can have several reasons. Unfortunately, the intrin- 
sic uncertainty of light curve modelling is hard to estimate 
and there is always some degeneracy between ^®Ni and to- 
tal mass, especially when radiation transport is treated in 
rough approximation. 

However, a s sumin g that the ^^Ni and total mass derived 
by ISauer et al.l l|2006l ) are correct, the discrepancy of the 
main properties can be explained by the over-simplified in- 
put model that was used in the calculation. It is well known, 
that mixing of the ejecta does have an important effect on 
the resulting spectra (also see above). 

While for the code comparison of RTM and NERO a 
one-zone model has been used, a real SN is certainly more 
complex. Apart from large scale asymmetries (which are not 
necessarily expected in SN 19941), the ejecta can be struc- 
tured on much smaller scales. If ^®Ni and other elements are 
separated, the ratio of Fe and other element lines changes. 
This, in turn, influences the estimate of the ejecta proper- 
ties. 

Mixing can influence the spectrum in two ways. First, 
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Figure 2. Synthetic spectra of tlie IWooslev et al. I | |1994) 13C 
model at 200 days after tlie explosion. The RTJ calculation is 
shown in black, while the red curve was produced using NERO. 
The agreement is excellent. 



2000 4000 6000 8000 10000 
Wavelength [A] 

Figure 3. Synthetic spectra of the IWooslev et al. I lUigl) 13C 
model at 400 days after the explosion. The RTJ calculation is 
shown in black, while the red curve was produced using NERO. 
There is some disagreement, especially in the Fe dominated region 
around 5000 A but in general the agreement is excellent. 



separation of Ni and other elements reduces the 7 heating 
in the non-radioactive zones. However, in the intermediate 
nebular phase this effect is weak, since the 7 opacity is low 
and positrons do not dominate yet. More importantly, when 
separated, carbon or oxygen rich zones cannot cool via Fe 
emission lines. This means, that separating ^''Ni from other 
elements in the nebular phase can lead to stronger emission 
of those elements than with perfect mixing, which may seem 
counter-intuitive on a first glance. 

To demonstrate this, we model SN 19941 using NERO 
again, this time separating Fe and O in several thin shells 
to simulate a separation of the ejecta on small scales. The 
derived ®^Ni and total mass changes strongly (see Table [TJ 
model 'C') and becomes more consistent with the hght curve 
estimate of lSauer et al.l l|2006l ). Of course, such an approach 
is highly degenerate and a broad variety of modelling results 
is possible. 

This causes an unfortunate situation. On the one hand, 
detailed knowledge of the mixing of the ejecta on large and 
on small scales is necessary to derive the main ejecta prop- 
erties. On the other hand, this information is poorly con- 
strained from observations and explosion models, especially 
on the small scales. Therefore, mixing poses a problem for 
deriving ejecta properties of stripped-envelope core-collapse 
SNe from modelling. Within these uncertainties, RTM seems 
appropriate to derive the main properties of SNe. For NERO 
calculations, more elaborate input models seem to be neces- 
sary to become consistent with the hght curve modelling. 



3.2 Comparison to observations 

To further test the reliabi lity of NERO, we calcu late synthetic 
spectra for a W7 model (iNomoto et al.|[l984 ). which is ex- 
pected to reproduce the spectra of 'normal' SNe la, although 
this has never been tested at intermediate e pochs. We com- 
pare our synthetic spec tra to SN la 2005cf (jCaravini et al.l 
I2OO7I : IWang et al.l|2009| ) , which can be regarded as a proto- 
type of 'normal' SNe la. The calculations are performed at 
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Figure 4. Synthetic spectra of the W7 model at 94 days after 
the explosion. The RTJ calculation is shown in black, while the 
red curve was produced using NERO. The agreement is reasonable, 
apart from a Fe III feature around 4600 A and a Co III feature 
around 5900 A. 



47, 94 and 338 days after the explosion (see Figures [S] |9] & 

At 47 days after the explosion the agreement between 
the synthetic and the observed spectrum is acceptable, given 
that Co coUisional data and forbidden lines are poorly 
known. It is interesting to note that the prominent feature 
at ~ 8500 A, which is usually thought to be a Ca 11 P-Cygni 
profile, could be strongly influenced by Co 11 emission in the 
observed SN la at this epoch. The coUisional data for Co 
II are very poor and it is not unlikely that we underesti- 
mate (or overestimate) the Co emission in our calculation. 
Since important Co coUisional data is missing, this remains 
a speculation. 

At 94 days after the explosion there is serious disagree- 
ment between the observed and synthetic NERO spectrum 
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Figure 5. Synthetic spectra of the W7 model at 338 days after 
the explosion. The RTJ calculation is shown in black, while the 
red curve was produced using NERO. There is some disagreement, 
but in general the agreement is good. 
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Figure 6. Synthetic spectrum of the W7 model at 338 days after 
the explosion. The RTM calculation is shown in black, while the 
red curve was produced using NERO. In general the agreement is 
good. 



around 4600 A. Interestingly, in the synthetic spectrum this 
feature is dominated by Fe ill (see Figure [9]), while another 
feature, which matches the observations well, at ~ 5900 A 
is dominated by Co ill. The RTJ spectrum of W7 (see Fig- 
ure [10]) fits the 4600 A much better, but under-produces the 
flux around 5900 A. Therefore, there seems to be a problem 
with the ratio of the Fe ill and Co ill line emission in both 
calculations. Although, one should not expect that the W7 
model can reproduce the spectra of SN 2005cf in all details, 
this could mean that our atomic data for Fe are inaccurate 
or that our approximations for example for the excited state 
ionisation cross-sections are too simple. In any case, the lack 
of reliable Co ll & ill data poses a serious problem for all SN 
la spectral calculations between maximum light and ~ 150 
days after the explosion before most of the Co has decayed 
into Fe. It has to be hoped that these data will be available 
in the near future. 



Figure 7. The spectrum of SN 19941 at about 159 days after 
the explosion is s hown in black. The spectrum is taken from 
ISauer et al.l l|2006ll and an identification of the most important 
emission lines can be found there. The RTM model is shown in 
blue, while the NERO models 'B' and 'C are shown in green and 
red, respectively. All spectra have been scaled by 2.2 ■ 10^^ cm'^ 
s ergs~^. All models seem to fit the observed spectrum equally 
well. Their properties are listed in Table [T] 



At 338 days after the explosion the flux is dominated 
by Fe emission lines and the agreement of the synthetic and 
the observed spectrum is good. At those epochs the nebular 
spectra of 'normal' SNe la are dominated by three prominent 
Fe features at roughly 4400, 4700 and 5300 A. The Fe 'tri- 
dent' is shaped by the ionisation state of the Fe core, which 
strongly depends on both the density of the core and the 
ratio of radioactive and stable iron. Also mixing with light 
and intermediate mass elements can influence the relative 
abundance of Fe ions. 

While the 4700 A feature is dominated by [Fe ill] emis- 
sion, the 4400 A feature is made from [Fe l] and [Fe ll]. In 
our synthetic spectrum this feature is underestimated and it 
may well be that it contains more contribution from [Fe l] in 
the observed SNe la than predicted in our simulation. It is 
important to note that very small fractions of Fe I (~ 0.01%) 
are sufficient to cause observable [Fe l] lines. Such small frac- 
tions of Fe I can survive even when Fe ill is present, strongly 
depending on photo-ionisation, recombination and possibly 
on charge exchange processes. This makes an accurate pre- 
diction of [Fe l] features difficult, at least at these epochs. 
The 5300 A feature contains both [Fe ll] and [Fe ill] and 
shows also a [Fe l] contribution. There is almost no Co emis- 
sion, except weak [Co ill] lines observed around 6000 A and 
[Co ll] emission at about 10000 A. Since we have no collision 
strengths for Co lines from the literature, their strength may 
be underestimated in our simulation. The RTJ and the RTM 
spectra of W7 at 338 days after the explosion are similar to 
the NERO spectrum and are shown in the previous section. 



4 DISCUSSION 

A comparison t o th e radiation transport code of 
IJerkstrand et all l|201ll ) has shown excellent agreement. 
This indicates that both codes work properly within the 
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Figure 8. The spectrum of SN 2005cf at about 47 days after the 
explosion l lWane et al.ll2009l 'l is shown in black (upper panel) or 
light grey (lower panel). The spectrum was scaled by a constant. 
The coloured curves were produced using NERO on W7. Upper 
panel: the red line shows the total synthetic flux. In general the 
agreement is good. It is important to note that the early-time 
ejecta are dominated by Co II & in, which have poorly known 
coUisional data. Lower panel: the flux of Co III (red), Co II (or- 
ange), Fe III (dark blue), Fe II (light blue) and Fe I (green) is 
shown separately. Most features are reproduced well. Disagree- 
ment is found around 4600 A (Fe II, Co ll), ~ 5900 A (Co III) 
and around ~ 8700 A, which is possibly Ca II and [Co ll] in the 
observed spectrum. Our synthetic spectra contain contributions 
from elements other than Fe or Co, which are not shown in the 
lower plot. 

uncertainties of the atomic data. This is an interesting re- 
sult, since they have been developed completely independent 
from each other and rely on a different numer ical approach. 

The comparison to the nebular code of iMazzali et al.l 

has also shown acceptable agreement, especially for 
pure Fe cores, which is important for SNe la. 

For SNe Ic, where several elements like C, O, Na, Mg, 
Si, S, Ca and Fe are important for the formation of the neb- 
ular spectra, we notice some differences. Most importantly, 
photo-ionisation influences the mass estimated for certain 
elements like Na and Mg. Also, there is some disagreement 
regarding the main properties of the SN core (total and ^^Ni 
mass) . 

It is well known, that mixing or a separation of the 



Figure 9. The spectrum of SN 2005cf at about 94 days after the 
explosion llGaravini et al.ll2007l ) is shown in black (upper panel) 
or light grey (lower panel). The spectrum was scaled by a con- 
stant. The coloured curves were produced using NERO on W7. 
Upper panel: the red line shows the total synthetic flux. In gen- 
eral the agreement is excellent apart from the region around 4600 
A. Lower panel: the flux of Co III (red), Co II (orange), Fe III (dark 
blue), Fe II (light blue) and Fe I (green) is shown separately. The 
synthetic flux exceeds the observed one at ~ 4600 A by a factor 
of ~ 3 and is dominated by [Fe III] emission. At the same time, 
the [Co III] feature at ~ 5900 A matches the observed spectrum 
well. 



elements on small or large scales in SN ejecta can have strong 
influence on the resulting nebular spectra. Since this paper 
intends to compare codes, we do not study this effect in 
detail. However, it was demonstrated that the uncertainty 
caused by the mixing of the ejecta is comparable to the 
uncertainty caused by using the different codes, at least for 
modelling SNe Ib/c. A treatme nt of the mixing problem i n 
SNe II has been prese nted by iKozma fc FranssonI (|l998h : 
Ijerkstrand et all (|201ll ). In 'normal' SNe la this problem 
is less severe, since the core is dominated by ^®Ni decay- 
products. 

By comparing synthetic la spectra to observations of the 
proto-typical 'normal' SN la 2005cf we have shown that the 
synthetic spectra produced with NERO look reasonable and 
are likely reliable within the uncertainties of the atomic data. 
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Figure 10. The s pectrum of SN 2005cf at about 94 days after 
the explosion (Gar avini et al. is siiown in black. The spec- 

trum was scaled by a constant. The red curve was produced using 
RTJ on W7. The reproduction of the 4600 A feature in the RTJ 
calculation is much better than in the NERO calculation. Other 
features, for example around 4000 A and 5900 A are reproduced 
better using NERO. In general, the agreement with the observation 
is excellent. 

Most of the observed spectral features have been identified 
to result from either Fe or Co emission. 

At epochs between 50 and 150 days after the explo- 
sion poorly known Co data pose severe problems for spectral 
modelling of SNe la. Atomic data are essential for calculat- 
ing SN spectra. Especially for electron collisions of all kinds 
the available data are often inaccurate or incomplete. 

Apart from time-dependent effects, NERO treats all the 
radiation transport effects commonly thought to be impor- 
tant for the formation of SN spectra in full NLTE. Therefore, 
NERO calculations are especially interesting for intermediate 
epochs, since so far SN spectral calculations at 50 — 200 
days after the explosion have been extremely rare. Also, the 
nebular phase between roughly 200 and 500 days after the 
explosion can be studied. 

Possibly, NERO could be used to calculate (pre-) maxi- 
mum spectra by imposing an estimated lower boundary flux 
at appropriate radii, a s it is done in photospheric codes (e.g. 
iMazzah fc Lucvlll993l ). With respect to purely photospheric 
codes, a treatment with NERO would include the effect of net 
emission above the lower boundary, which could then be set 
at lower velocities than in previous approaches. This may 
increase the accuracy of the (quasi) photospheric approach 
considerably. We plan to investigate this possibility in the 
near future. 



5 SUMMARY AND CONCLUSION 

In this paper we presented a new NLTE radiation transport 
code, which can be used to calculate synthetic spectra for all 
types of SNe at intermediate and late epochs. Our treatment 
of intermediate epochs opens a new window for SN spectral 
analysis. Currently, NERO is working in spherical symmetry, 
but a three-dimensional version may be available in the fu- 
ture. In its one-dimensional version the code can be used for 
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Figure 11. The s pectrum of SN 20 05cf at about 338 days af- 
ter the explosion l lWane et al.ll2009l ) is shown in black (upper 
panel) or light grey (lower panel). The spectrum was scaled by 
a constant. The coloured curves were produced using NERO on 
W7. Upper panel: the red line shows the total flux. In general 
the agreement is good. Lower panel: the flux of Co in (red), Co 
II (orange), Fe III (dark blue), Fe II (light blue) and Fe I (green) 
is shown separately. Co features can be observed around 6000 A 
and around 10000 A. Their strength may be underestimated in 
our calculation. 

spectral modelling of observed SN spectra or for calculating 
synthetic spectra of (approximately) spherically symmetric 
SN explosion models. 
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